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THE MID-INFRARED TULLY-FISHER RELATION: SPITZER SURFACE PHOTOMETRY 
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ABSTRACT 

The availability of photometric imaging of several thousand galaxies with the Spitzer Space Telescope 
enables a mid-infrared calibration of the correlation between luminosity and rotation in spiral galaxies. 
The most important advantage of the new calibration in the 3.6/im band. IRAC ch.l, is photometric 
consistency across the entire sky. Additional advantages arc minimal obscuration, observations of flux 
dominated by old stars, and sensitivity to low surface brightness levels due to favorable backgrounds. 
Through Spitzer cycle 7 roughly 3000 galaxies had been observed and images of these are available 
at the Spitzer archive. In cycle 8 a program called Cosmic Flows with Spitzer has been initiated 
that will increase by 1274 the available sample of spiral galaxies with inclinations greater than 45 
degrees from face-on suitable for distance measurements. This paper describes procedures based on 
the photometry package Archangel that are being employed to analyze both the archival and the new 
data in a uniform way. We give results for 235 galaxies, our calibrator sample for the Tully-Fisher 
relation. Galaxy magnitudes are determined with uncertainties held below 0.05 mag for normal spiral 
systems. A subsequent paper will describe the calibration of the [3.6] luminosity— rotation relation. 
Subject headings: infrared: galaxies; galaxies: photometry, distances 



1. INTRODUCTION 

The Cos mic Microwave Back ground (CMB) tempera- 
ture dipole ([Fixsen et al .11 19961 ) is usually interpreted as 
a motion of our Galaxy of over 600 km/s but the con- 
siderable majority of the posited motion is developed on 
large scales with origins that are still poorly understood. 
Our overarching goal is to measure distances, hence 
parse departures from the mean Hubble expansion, on 
scales extending to 200 Mpc. We are gathering distance 
measures from a multitude of methods and contributors 
within a program that we arc calling Cosmic Flows 
(see Appendix). Of particular importance for us are 
distances accrued from the correlation between the 
rotation rate of a gala xy and its luminosity , the Tully- 
Fisher Relation: TFR (jTullv k Fisherll 19771) . There are 
methodologies that provide distance estimates that are 
individually more accurate but an abiding advantage of 
the TFR is applicability to a large fraction of all galaxies 
over a wide range of environments and distances. There 
is the prospect of utilizing the TFR to obtain distances 
to several tens of thousands of galaxies out to redshift 
z ~ 0.05 (~ 200 Mpc) . 

As steps toward the accumulation of an appropriate 
data set of distances, our Cosmic Flows Large Program 
on the 100m Green Bank Telescope and complemen- 
tary southern observations on the Parkes Telescope 
(|Courtois et al.ll2009l l2011bO are providing us with qual- 
ity rotation information from HI line profiles, and the 
merging of our own a nd literature optical photometry 
(jCourtois et al.ll2011a! ) provides the other element that 
has permitted a modern re-calibration of the TFR at 
an optical band (|Tullv k Courtois! [2012]) . The data 
accumulated in support of this program, including the 
new material to be discussed in this paper, are made 
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available at EDD, t he Extragalactic Distanc e Database, 
acces sed online at http://edd.ifa.hawaii.edu (jTullv et al.1 
I2009T ). 

It has long been appreciated that photometry in the 
infrared may offer advantages because of reduced extinc- 
tion and because infrared flux arises in large measure 
from old stars that should dominate the inventory of 
baryonic mass (jAaronson et al.1 11979D . Photometry at 
K s band from the Two Micron All S ky Survey (2M ASS) 
has been used with the TFR (jKarachentsev et all 
|2002[ ). However a major concern with observations 
in the infrared from the ground is high and variable 
sky foreground. Much of the flux from galaxies lies in 
extended components with surface brightnesses that 
are well below the ground based sky level. Flux at 
the extremities of galaxies is lost and very low surface 
brightness galaxies are not even seen. 

Observations from space removes the problem of the 
high contamination by the Earth's atmosphere. We 
have initiated a sub-program that we call Cosmic Flows 
with Spitzer (CFS) w ith NASA's Spitzer Space Telescope 
(| Werner et al.ll2004l ). Observations have begun in cycle 
8 during the post-cryogenic period to obtain wide-field 
images of galaxies with IRAC, the InfraRed Array Cam- 
era, in Channel 1 (ch.l). The present paper describes 
our reduction and photometry analysis procedures. We 
will discuss the transformation steps from raw Spitzer 
Post Basic Calibrated Data (PBCD) obtained from the 
Spitzer archive to the parameter needed for the TFR cal- 
ibration: apparent [3.6] magnitudes. The photometry is 
carried out with a Spitzer-adapted version of Archangel 
(jSchombertJ [20071 ISchombert fc Smith! [2011 . We will 
discuss the corrections to be made to apparent magni- 
tudes and conclude with a discussion of uncertainties. A 
subsequent paper will discuss the calibration of the TFR 
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at 3.6/im. 



2. [3.6] BAND DATA 



Once cryogens were exhausted on Spitzer Space 
Telescope useful observations were restricted to IRAC 
ch.l (centered at 3.55/im which we round to 3.6/im) 
and eh. 2 (4.5/mi). For the purposes of measuring 
distances with the TFR the two channels are highly 
redundant. Given a choice with the availability of 
finite observing resources between more galaxies in 
one band versus fewer galaxies in two bands we chose 
more galaxies in one band. In our CFS program we 
concentrate on IRAC ch.l observations in the 3.6/im 
window that give us magnitudes [3.6] in the AB system. 
This wind ow provides observa tions with minimal dust 
extinction ([Draine fc Ledfl984l) . It lies at a minimum o f 
the zodiacal background radiation (|Ootsubo et al.llT998D . 

Figure [T] provides examples of th e spectral energ y 
distribution (SED) of spiral galaxies (jSilva et al.l 119981 ). 
The Spitzer [3.6] band lies on the Rayleigh- Jeans tail 
of the SED of normal populations of stars, not yet 
strongly affected by flux from warm dust that starts to 
become a factor at longer wavelengths than 4/im. The 
discrete spectral features seen in the SED arise from 

Polycycl ic Aromatic Hydrocarbon (PAH) molecules 

(|Tielensll2008h . The highest frequenc y PAH, at 3 3/jm is 
contained within the [3.6] bandpass. iMcidt et al.l (|2012[ ) 
have investigated the impact of various contributors to 
flux in the 3.6/im window with 6 representative spiral 
galaxies observed with the Spitzer S4G program. They 
find contributions from hot dust and PAHs together 
contribute 9 ± 4% of the global flux in the 3.6/im 
band, intermediate age asymptotic giant branch and red 
supergiant branch stars contribute 3 ± 2% of the global 
flux, and the rest, the great majority, is contributed 
by old stars, predominantly K and M giants. These 
non-stellar and young stellar contributions should only 
slightly degrade the correlation between old stars and 
mass in normal spirals. 
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Fig. 1. — Comparative SED for spirals of types Sa (red, dashed), 
Sb (black, solid), and Sc (blue, dotted). The Spitzer [3.6] passband 
is illustrated along with the wavelengths associated with B,R,I,K 
bands. The relative scales of the SED are offset to match at 0.8/im. 

Using Spitzer IRAC ch.l, a point spread function with 
mean FWHM 1.66" is sampled with 1.2" pixels. The 



field of view is 5.2', adequate to encompass most galaxies 
to beyond twice 0I25, the diameter at a B isophote of 
25 mag per square arcsec. Larger galaxies require 
mosaics. Integrations with the CFS program involve 
the combination of 8 x 30 second slightly dithered 
exposures for a total of 4 minutes per field. As will be 
discussed, these integrations provide images that probe 
somewhat fainter limits than most ground-based optical 
photometry programs and much fainter limits than 
ground-based infrared photometry programs. Spitzer 
surface brightness levels reach ten magnitudes below 
typical ground-based infrared sky levels. No existing 
near-infrared ground survey achieves the accuracy 
obtained with Spitzer Space Telescope. The outstanding 
advantages of space observ ations are backgro und stabil- 
ity and all-sky consistency (jFazio et al.ll2004T ). 

Our cycle 8 post-cryogenic program CFS avoids rep- 
etition of earlier Spitzer observations. Archival infor- 
mation is used where available. Major contributions 
from earlier programs come from SINGS, the Spitzer In- 
frared Nearby Galaxies Survey (|Dale et all [ 20051 [2001 
and L VL, the Local Volume Legacy survey ([Dale et al.l 
120091 ) carried out during the cryogenic phase, then 
S 4 G, the Spitzer S urvey of Stellar Structure in Galaxies 
(ISheth et al.H2010h. and CHP , the Carnegie Hubble Pro- 
gram ( Freedman et al.ll2011[ ). subsequently carried out 
during the post cryogenic phase. Smaller programs un- 
dertaken during the cryogenic phase supply us with a few 
more fields. These data are available forpublic use at the 
Spitzer Heritage Archive website (SHAQ). The variety of 
source programs introduces variations in the details of 
the acquisition, particularly affecting total integrations, 
dithering procedures, and the extent of fields referenced 
to d,25- However, with all the data that will be consid- 
ered the fields are large enough and the exposure times 
are long enough that at most only a few percent of the 
light from a target is lost. 

3. PHOTOMETRY 

3.1. Surface Photometry 

Large numbers of pixels complicate simple parameter 
extractions. A galaxy is spread over a large area of the 
sky. At some point outer pixels have more sky lumi- 
nosity (zodiacal light and background contaminants) 
than galaxy luminosity. Setting the "sky" dominates 
the total magnitude error budget. An analysis of a 
large galaxy (extending across many pixels) requires 
surface photometry involving fits of isophotes, lines 
of constant luminosity. Isophotes are often set to 
be ellipses (jMilvang- Jensen fc Jorgense 3 H999T) . Our 
interest is with spiral galaxies with types typically 
between Sa and Scd. A well behaved spiral is approx- 
imated by an oblate spheroid that appears circular 
when viewed face-on and projects to an ellipse when 
viewed toward edge-on. Galaxy 2D images described 
by elliptical isophotes can be summed in annuli to 
reduce to a ID description. Then, the ID profiles are 
fitted by various functions in order to extract the radial 
surface brightness (SB) distribution, global structure or 
geometrical characteristics, spatial orientation, stellar 

1 http: / /irsa. ipac.caltcch.edu/data/SPITZER/docs/spitzcrdataarchivcs/ 
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populations, charact eristics of dust, etc. To obtain 
apparent magnitudes, Ide Vaucouleursl (I1977D introduced 
the growth curve, a plot of magnitude within a radius as 
a function of radius. With an adequate signal to noise 
ratio, it could be enough to place large apertures around 
galaxies and sum the total amount of light, minus 
the sky contribution. In practice, a galaxy luminosity 
distribution decreases towards larger radii so larger 
apertures catch more galaxy light but also introduce 
more sky noise. Some light is inevitably lost below 
the sky level. Isophotal intensities associated with the 
galaxy light at large radii are sensitive to the sky setting. 
Restriction to a smaller radius leads to underestimates 
of total light. The problem is that galaxies do not have 
discrete edges. 

It is never possible to measure 100% of the light of 
a galaxy. Measurements are made to an isophotal level 
dictated by telescope optics, detector, exposure times, 
and sky brightness. Different authors measure magni- 
tudes to different isophotal levels then often extrapolate 
to total magnitudes. Our interest is with spiral galaxies. 
These galaxies characteristically decay exponentially in 
luminosity with radius. In an ideal case, light contained 
within a specified isophotal level is a simple function 
of the disk central SB and the exponential decay scale 
length. To extrapolate in such a case one can assume 
that the light at large radii falls off like an exponential 
disk with a central SB and scale length characterized by 
a fit to the main body of the galaxy. The estimated con- 
tribution lost below the sky level can be adde d to what is 
obser ved to give an extrapolated magnitude (jTullv et al.l 
[19961) . 

3.2. Archangel 

iSchombertl (|2007l ) developed Archangel, a flexible 
tool for galaxy surface photometry built of a combi- 
nation of FORTRAN and Python routines. Archangel 
performs procedures such as: 1) masking of stars and 
flaws, 2) ellipse fitting at expanding radii from the 
galaxy center, 3) compression of 2D information into 
ID SB and magnitude growth curves as a function of 
radius, and 4) extrapolation via fits to the magnitude 
rowth curve at larg e radii involving rational functions 
Schombert k, Smith! I2012D . Position angles and el- 



lipticities are freely determined at each radial step in 
the development of the growth curve. At large radii 
noise dominates and position angle and ellipticity are 
frozen for the remaining outward steps in radius. The 
program provides flexibility in where these parameters 
become frozen and allows that they may be frozen 
at all radii. Total magnitudes, the most important 
product of this analysis, are found to be negligibly 
affected by position angle and ellipticity details at 
intermediate radii. Comparisons with alternative pho- 
tometry are discussed in section 5.1. See Figure [2] for 
an example of masking and ellipse fitting with Archangel. 

To obtain magnitudes in the AB system from the 
archival PBCD we use magnitude and VEGA/AB con- 
version parame ters from the IRA C Instrument hand- 
book and from lCaputi et akl ()2006l ) respectively. At op- 



tical bands it is common practice to quote magnitudes 
in the Vega photometric system but working in the mid- 
infrared it is more useful to use the AB system. In 
instances where comparisons are made between optical 
and mid-infrared, w e use the following transformations 
dFrei fc Gunnl[l99l : 

B{Vega) = B{AB) +0.163 

R c {Vega) = R C {AB) - 0.117 

I c {Vega) = I C (AB) - 0.342 

[3.6] (Vega) = [3.6] (AB) - 2.785 
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Fig. 2. — Output of the Archangel software showing the ratio 
b/a and the position angle in the two top panels and the fitted 
ellipses and the masking in the two bottom panels for PGC41729 
= NGC4522. 

A significant source of uncertainty arises from the 
setting of the sky level (see section 5). In Archangel 
the sky is taken as the median of sky boxes placed 
around the galaxy. This method gives rea listic initial 
sky background estimates (|Hall et al.l |20 111) . If targets 
are modest in size there is reasonable control of the 
sky level. If the sky is set properly then the magnitude 
growth curve should go asymptotically flat at large 
radii. One can also evaluate the sky setting by looking 
at the SB as a function of radius. SBs are not expected 
to flare or drop precipitously at the sky level, although 
such occurrences are no t phenomenologicall y excluded 
(jMacArthur et al.l l200l lErwin et al.l l200l . Visual 
inspections of the magnitude growth curve and SB 
dependence with radius ensure an optimal sky setting. 
Fortunately, sky values are low in Spitzcr data even if we 
will show in the last section that this problem remains 
our major source of uncertainty. 



http:/ /irsa. ipac.caltech.edu/data/SPITZER/docs/irac/ 
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An issue related to the sky problem is the matter of 
the terminal radius of an analysis. A limit to the fitting 
process can be imposed by signal-to-noisc considerations. 
Integration times permit us to reliably reach a radius 
026.5 at the isophotal level 26.5 mag arcsec~ 2 in the [3.6] 
band. We try to extend the ellipse fitting to 1.5<226.5 in 
the [3.6] band. A goal of the program is to assure that the 
ellipse fitting extends to at least 1.1 026.5, with mosaics 
if necessary. This [3.6] band dimension is not available 
before the observation so we rely on a substitution found 
to be comparable based on the B band diameter c?25, 
requiring that the observed area extend to a radius 1.5e?25 
(|Sheth et al.l[2010h . 

3.3. ID Fits and Parameter Extractions 

The mean SB in magnitudes per square arcsec in an 
annulus at radius depends on the mean flux in a pixel 
at that radius Ffa) and the mean sky flux in a pixel S: 



H(n) = -2.51og( 



0.6 ' 



21.585 



(1) 



where the constant in the denominator provides a 
conversion from pixels to arcsec. 
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Fig. 3. — Output of the Archangel software showing the SB profile 
for PGC41729=NGC4522. The red, green, and blue lines identify 
the radii enclosing 20%, 50%, and 80% of the total light respectively 
and the surface brightnesses at those radii. The dashed straight 
line illustrates an exponential fit to the disk between the radius 
enclosing 50% of the light and the isophotal radius (126.5- The 
insert provides a [3.6] image of the galaxy with an annulus at the 
isophotc 26.5 mag arcsec -2 . 



Archangel allows for a description of the run of SB 
with radius as the sum of disk and bulge components. 
Instead, we choose to restrict to disk fits only. With 
multiple component fits there are frequently trade-offs 
such that the overall fit may be satisfactory but the 



physical meanings of parameters are ambiguous. Usually 
the dominant radial SB characteristic of spiral galaxies 
is an exponential decay of projected luminosity with 
radius. Deviations are most frequently seen toward 
the center where a bulge may become dominant. It 
is beyond the scope of this program to dissect galaxy 
images into detailed morphological components because 
such dissection has negligible effect on the product that 
most interests us: total magnitudes. We restrict fitting 
to a rough characterization of the exponential fall-off. 

The radial run of SB in an exponential disk, /^(r), has 
a behavior described by an c-folding scale length, a, and 
central SB, 



t*d(r) = Mo + 1.0857(r/a). 



(2) 



As a procedure, we determine the disk parameters with 
a fit over the range of radii from the half-light radius 
a e to the isophotal limit radius 026.5- This range is 
modified if a visual inspection indicates the need. An 
example of a SB product is illustrated in Figure [3] 

An example of a magnitude growth curve as a function 
of semi-major axis is shown in Figure |4j The light from 
each succeeding annulus contributes to the (negatively) 
increasing magnitude with increasing radius. If the 
sky value is properly set then the growth curve will 
asymptotically flatten. Should the curve turn over it 
would be inferred that the sky level is set too high - 
flux from the galaxy is being attributed to the sky and 
being removed. Conversely, the sky set too low causes 
flux from the sky to be attributed to the galaxy and the 
growth curve will fail to flatten. 

The Spitzcr photometry is sufficiently deep that magni- 
tudes in the growth curve approach the total magnitude 
of the galaxy. One way to extend to the total magnitude 
uses the procedure built into Archangel based on interpo- 
lations and extrapolations with rational functions. Such 
functions have a wide range in shape and have better in- 
terpolating properties than polynomial functions. They 
suit data where an asymptotic behavior is expected. The 
quadratic form/quadratic form, meaning a degree of 2 in 
both numerator and denominator is the simplest choice. 
The asymptotic magnitude is Ci/ca where c\ and C2 are 
the second order coefficients of the numerator and de- 
nominator respectively. However, rational functions are 
non-linear. They can produce vertical asymptotes due 
to roots in the denominator that are to be ignored. Fit 
uncertainties are given by the standard error of the esti- 
mate, SEE: 



SEE = 



\ 



1 ™ 

-^2(m(ai) f it 



m(ai) 



measure^ 



(3) 



Given a growth curve as seen in Figure [4] it is straight- 
forward to define the useful parameters 020, a e , and a^o 
enclosing 20%, 50%, and 80% of the light respectively. 
The associated magnitudes and semi-major radii are 
illustrated in Figures 2—4. Other products are the 
average SB within a e and 020 and a concentration index 
C's2 = Q8o/o20- Table [1] gives the parameters that are 
extracted for the galaxy used as an example in Figs. 
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2—4 and illustrates what is seen in a single row in the 
catalog Spitzer [3.6] Band Photometry at EDD, the 
Extragalactic Distance Databasfl Capabilities within 
EDD allow a user to link to other catalogs, thereby 
accessing all manner of information about each target. 

We provide three magnitudes that approximate the 
global magnitude of the galaxy: [3.6]26.5, an isophotal 
magnitude that directly measures the light within a re- 
liably attained radius, [3.6]t {, a 'total' magnitude given 
by the asymptote of the Archangel rational function ex- 
trapolation, and [3.6] ea; t, an 'extrapolated' magnitude as- 
suming a continuation of the exponential disk beyond the 
radius of the isophote 26.5 mag arcsec -2 . The relative 
merits of these magnitudes will be discussed in a later 
section. 

PGC41729_2 



Mag: 1 1 .970 +/- 0.003 



10 



11 



! 12 



I 



13 T 



ft 
ft 
ft 



16 



17 



.ft* 



og(semi-major axis) (arcsec) 

Fig. 4. — Output of the Archangel software showing the growth 
curve for PGC41729=NGC4522. The incremental growth of the 
apparent magnitude of the galaxy with radius is shown by the 
progression of stars. The fit providing an extrapolation to a total 
magnitude is generated over the domain of the red stars. The level 
of the total magnitude is shown in magenta. The 80%, 50%, and 
20% light enclosure levels are shown in blue, green, and red. 
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TABLE 1 

Extracted photometry parameters. (1) Principal Galaxies 
Catalog number, (2) common name, (3) date of Spitzer 
observation, (4) nominal total integration, seconds 
(actual time collecting photons somewhat less), (5) (126. 5: 
major axis radius at isophote 26.5 mag arcsec -2 , (6) 

[3.6]26.5: AB MAGNITUDE WITHIN 026.5, (7) [3.6]tot: TOTAL AB 
MAGNITUDE FROM RATIONAL FUNCTION ASYMPTOTE, (8) <T m : RMS 
DEVIATIONS, RATIONAL FUNCTION FIT, (9) [3.6]eict: TOTAL AB 
MAGNITUDE BY EXTRAPOLATING FLUX BEYOND £126.5 ASSUMING 
CONTINUANCE OF EXPONENTIAL DISK, (10) M0: CENTRAL DISK 
SURFACE BRIGHTNESS FROM INWARD EXTRAPOLATION OF DISK FIT, 
MAG ARCSEC -2 , (11) a: EXPONENTIAL DISK SCALE LENGTH, 
ARCSEC, (12) b/a: RATIO OF MINOR TO MAJOR AXES, (13) cr b j a : 
UNCERTAINTY IN AXIAL RATIO, (14) PA: POSITION ANGLE OF 
MAJOR AXIS, DEG. (15) ago: MAJOR AXIS RADIUS OF ANNULUS 
ENCLOSING 80% OF TOTAL LIGHT, ARCSEC, (16) M80 ! SURFACE 

BRIGHTNESS AT ago , MAG ARCSEC -2 , (17) a e : 'EFFECTIVE 
RADIUS', MAJOR AXIS RADIUS OF ANNULUS ENCLOSING 50% OF 
TOTAL LIGHT, ARCSEC, (18) /i e : SURFACE BRIGHTNESS AT a e , MAG 
ARCSEC -2 , (19) < Me >: AVERAGE SURFACE BRIGHTNESS WITHIN 
a e , MAG ARCSEC -2 , (20) £120 : MAJOR AXIS RADIUS OF ANNULUS 
ENCLOSING 20% OF TOTAL LIGHT, ARCSEC, (21) /^20 : SURFACE 
BRIGHTNESS AT O20, MAG ARCSEC -2 , (22) < [120 >■ AVERAGE 
SURFACE BRIGHTNESS WITHIN £X20 , MAG ARCSEC -2 , (23) C7g2: 
CONCENTRATION INDEX, ago/a20, (24) SPITZER PROGRAM LINK. 



4. CORRECTIONS 
Corrected apparent magnitudes [3.6] b,fc,I ' a for Spitzer 



IRAC ch.l data are given by: 



[3.6] 



b,i,k,a 



[3.6] -A 



[3.6] 



A [3.6] _ A [3.6] + A [ 3 .6] (4) 



with apparent magnitude [3.6] output from Archangel, 
galactic extinction correction Aj 3 ' 6 ', internal extinction 
correction A? 3 ' 6 ', /c-correction j4.J. 3 ' 6 ', and aperture correc- 
tion ^L 3 ' 6 '. We shall describe these terms in the following 
sub-sections. 

4.1. Galactic extinction correction 

Galactic extinction depends only on object coor- 
dinates and observational wavelengths. The In- 
fraRed Science Archive (IRS A) provides an online 



tool at http://irsa.ipac.ca ltech.edu/ap plications/DUST/ 
with 100/iin cirrus maps (jSchlegel et al WlMi) that sup- 
plies us with differential redde nings, E(B — V). We use 
the correction term given by iCardelli et alj (I1989T ) ac- 
counting for a small shift to the centroid of the Spitzer 
passband: 



A 



[3.6] 



R [3S ]E(B - V) 



(5) 



with i?[3.6] = 0.20. Galactic extinction magnitude 
corrections at [3.6] are only 9% compared to at Iq and 
4% of the corrections at B. Corrections at latitudes 
above 15° are almost always 0.05 mag or less, with 
uncertainties ~ 0.01 mag. 



12 


13 


14 


15 


16 


17 


18 


19 


b/a 




PA 


£»80 


M80 


a e 


Me 


< Me > 


0.26 


0.01 


34 
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36 
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3 http://edd.ifa.hawaii.edu 



4.2. Internal extinction correction 

Internal extinction is usually the greatest concern. For- 
tunately, in the infrared suc h extinction is very small. 
iGiovanelli et all (|1995l 119971 ) showed that there" is a lu- 
minosity dependence to galaxy internal obscurations. 
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iTullv et al.l (|1998ft confirmed and provided an alternative 
description of the effect. There is a subtle problem be- 
cause absolute magnitudes ar e not known a prio ri. They 
are a product of the analysis. ITullv et al.1 ()1998f ) framed 
magnitude corrections in term of a distance-independent 
surrogate, the line width parameter, W^x- Accordingly, 
the internal extinction correction can be written: 

,[3.6] 



7[3.6]log(a/6) 



where 



7[3.6] 



7[3.6] =0.10 + 0.19(log - 2.5) 



(6) 



(7) 



if W^x > 94 km s _1 and 7[ 3 6 ] = otherwise. W^ x 
is a measure of twice the maximum rotation rate of a 
galaxy derived from W m 5o, the HI profile width at 50% 
of the mean flux within th e velocity range encompassin g 
90% of the total HI flux (jCourtois et al.l [2001 I2011bl ). 
The measure includes a deprojection appropriate for the 
inclination i (see the above references for the derivation 
of inclinations) . 

There is an advantage to this formulation of the 
internal extinction. If the inclination is underesti- 
mated, log(a/&) is underestimated driving A^ 3 ' 6 ' low 



but then W l mx 
hence A' 3 ' 6 ' up 



is overestimated which drives 7[3.6], 

The two terms in A[ 3 ' 6 ' are affected 
in opposite directions. Regardless, internal absorption 
corrections are always small, rarely reaching 0.1 mag. 
Uncertainties in these corrections are less than 0.02 mag. 



4.3. K-correction 

The fc-correction (IQke k, Sandag c 1968) is small at the 
redshifts we encounter. IHuang et al.l ([2007D show a linear 
dependence of the fc-correction with redshift at 3.6/im. 
This linear dependence is independent of the galaxy type 
at small redshifts at this position on the Rayleigh- Jeans 
tail of the spectral energy di stribution of s tar lig ht. We 
use the low-z formulation by IHuang et al.l (|2007D : 



A 



[3.6] 



-2.27 z 



(S) 



with z the galaxy redshift. Uncertainties are at the level 
of 0.01 mag or less. 



4.4. Aperture correction 

The fourth and last adjustment is the aperture cor- 
rection. Aperture corrections are required for extended 
source photometry with Spitzer (e.g. galaxies) because 
their absolute calibrations are tied to point sources with 
IRAC observations. There is extended emission from the 
Point Spread Function outer wings, and the scattering 
of the diffuse emission across the focal plane that is 
captured by the extended source photometry but not by 
the calibrations on point sources. Since the photometry 
is normalized to 12" apertur es a correct io n mu st be 
applied for large apertures (jReach et al.l (|2005l ) and 
IRAC Instrument Handbook). The following correction 
is recommended^. For an effective aperture radius r in 

4 http:/ /irsa. ipac.caltech.edu/data/SPITZER/docs/irac/ 

iracinstrumenthandbook / 



arcseconds, the ch.l IRAC extended source aperture 

Correction is f IRAC true = flRAC measured X (Ae~ T + C) 

where A = 0.82, B = 0.37 and C = 0.91. The extended 
source aperture correction in magnitudes is: 



A [3.6] = _ 2 .51og(Ae- rB +C). 



(9) 



The average correction for galaxies of interest to our pro- 
gram is 0.10. The variations on this correction from 
source to source for our galaxies, which are typically 
larger than 1', is 0.01 mag and 10% relative uncertainties 
in the adjustment are negligible. 

5. UNCERTAINTIES 

An extremely important virtue of the Spitzer [3.6] 
band photometry is the robustness of the luminosity 
measurements (a) with uniformity across the sky, (b) 
with inclusiveness of target light because of the sensitiv- 
ity, and (c) because adjustments are small. There was a 
discussion of uncertainties associated with the adjust- 
ments in the last section and it can be summarized that 
as long as sources are not in extremely obscured regions 
of our Galaxy (At < 1) then the global uncertainty in 
adjustments is at the level of 0.03 mag or less, with 
internal absorption within sources dominant in the error 
budget. The IRAC handbook gives a 2 — 3% error on 
the absolute flux calibration (excluding the aperture 
correction), but more importantly for this program, 
claims photometry is repeatable across the sky at the 
1% level. 

Among our parameters we determine isophotal, 'total', 
and 'extrapolated' magnitudes. The latter two both 
approximate the global magnitude, the 'total' from the 
rational function asymptote of the growth curve and 
the 'extrapolated' from the extension of the exponential 
disk fit beyond the radius of the isophotal magnitude. 
By construction, [3.6]26.5 is fainter than [3.6] ex t and 
should be fainter than [3.6]t D t. The average difference 
< [3-6]26.5 — [3.6] e2 - t >= 0.016 mag corresponds to a 
typical disk fit of 6.2 exponential scalelengths at the 26.5 
mag/arcsec~ 2 isophote. The typical uncertainty in this 
extrapolation is below 0.01 mag except if the target is 
extremely low surface brightness. SB profiles of spirals 
can depart from a pure exponential at large radii, either 
with flares or truncations a nd because of the interplay 
between bulges and disks (jKentl I1985D . Yet because 
such a large fraction of the flux is captured by the deep 
Spitzer integrations the differences between measured 
isophotal and extrapolated magnitudes are so small as 
to leave little room for uncertainty in the extrapolation. 

By comparison, < [3.6]26.5 — [3.6] tot >= 0.007 mag, 
that is, [3.6] tot is fainter than [3.6] ea;4 by 0.009 mag on 
average. The rms scatter is 0.018 mag between these 
alternative measures. The differences are primarily 
due to a slight instability in the rational function fits. 
We give preference to the exponential disk extrapolations. 

We turn to what is probably the largest source of 
error, the setting of the 'sky' level. With observations 
in space at [3.6] band this noise level is dominated by 
diffuse zodiacal light and discrete high redshift galaxies. 
The discrete contaminants can be easily seen to very 
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faint levels in regions beyond the galaxy. They are less 
easy to see and exclude if they are superimposed on the 
target galaxy. A major task before running a surface 
photometry analysis is the removal of contaminants 
like foreground stars and background galaxies. Our 
approach is to not be too aggressive with the removal of 
contaminants. We remove contaminants as best we can 
on the target and remove contaminants in the adjacent 
sky to the same level, leaving in place fainter sources 
since such sources must also be hidden within the galaxy. 

It was described in the section on Archangel pho- 
tometry that sky settings were established from the 
median of pixel fluxes in boxes placed around the galaxy 
and validated by the nature of the magnitude growth 
curve (it should go asymptotically flat) and the surface 
brightness profiles (flares or cutoffs as noise dominance 
is approached arc suspicious but not considered a 
conclusive sign of bad sky setting) . In order to generate 
a quantitative test of the effects of sky variance we have 
run Spitzer-adapted Archan gel on 235 galax ie s, ou r 
calibrator sample defined in iTullv fc Courtoisl <|2012f ). 



Our calibrator sample contains only a few LSB and 
irregular galaxies. 

A first run gives us the sky value 5*0 and its uncer- 
tainty. We run Archangel two more times with sky values 
of So ± cr s k y respectively for each one of our selected 
galaxies. This gives us three extrapolated magnitudes 
that we call [3.6]o, [3.6] + , and [3.6] _. Figures [5] and |6] 
show the variation of (| [3.6]n — [3.6]_ | + 1 [3.6]o — [3.6] + |)/2 
as functions of type and apparent magnitude. These 
plots show the sensitivity to the choice of sky value and 
that this sensitivity becomes particularly acute for low 
surface brightness systems. The galaxies of type Sd, 
identified in the plots as low surface brightness galaxies 
(LSB), and the magellanic irregular galaxies are clear 
outliers. There are also 3 galaxies with very bright 
objects nearby that could influence the sky level with 
stray light and might explain their position in Figure [6l 
These 3 are retained in the calculations of the offset and 
scatter since they are typical spiral galaxies but the LSB 
and irregular galaxies are excluded. Magnitude uncer- 
tainties due to the sky error are of the order 0.04 ± 0.02. 
Uncertainties with LSB and irregular galaxy magnitudes 
tend to be more important which is understandable. 
LSB have surface brightness values closer to that of the 
sky and irregular galaxies not only tend to be LSB but 
in addition might not be well described by elliptical 
isophotcs. Changing the sky value a little might change 
the measured flux considerably toward the external 
part of such galaxies. LSB and irregular galaxies apart, 
Figure [6] indicates that uncertainties do not strongly 
increase at fainter magnitudes. 

Next we test for uncertainties in magnitude due to the 
sky against the sky value itself, as well as against the 
isophotal semi-major axis in the [3.6] band 026.5, incli- 
nation from face-on, and apparent area defined as the 
area of the ellipse at 026.5 to see if any trends exist. The 
results in Figures [7] and [5] show no correlation. In Figure 
El we can see that the uncertainty in magnitude does not 
depend on the sky value. We checked for a dependence 
on sky uncertainty and find no correlation. These re- 



CO 



0.3 



0.2 



0.1 



v 0.0 



co 



-0.1 



* Object close 
RMS = 0.02 
Offset = 0.04 
No.: 203 



* > 
+ i 

+ + + : 

+ + i i i = 


s + 

I i 
t * + 












+ + r 1 1 








i t + + 










+ 



S0°S0"S0/aSa Sab Sb Sbc Sc Scd SdSdmSm Im 
Morphological Type 
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tively excludes types Sd and later. Three cases with contamination 
from nearby bright objects are indicated by asterisks. The scatter 
is asymmetric about the mean since an absolute value difference 
from the fiducial value cannot be less than zero. 
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Fig. 6. — Variation of magnitude uncertainty as a function of 
magnitude. Type Sd systems, here referred to as low surface bright- 
ness galaxies, are represented by squares and types Sdm-Sm-Im 
irregular galaxies are represented by triangles. Asterisks locate 
galaxies with a very bright object close to them. The mean offset 
and scatter lines have the same meaning as in the previous figure. 

suits suggest that the total and extrapolated disk appar- 
ent magnitudes are adequate (we do not show the plots 
for both magnitudes here as they are very similar). In 
any case, the highest sky values are relatively moderate 
(< 0.20 MJy sr _1 ). One can also notice that sky values 
and sky uncertainties are not correlated, evidently a re- 
flection of the relative uniformity of background across 
dimensions of 5 — 10 arcmin. Structure in the background 
could be a worse problem when the sky setting is very 
low. Perhaps it is a surprise that the uncertainty is not 
proportional on the galaxy apparent area (Figure [51 bot- 
tom) . The more pixels that are affected by setting a new 
sky, the more the magnitude might change. In any case, 
these tests indicate that magnitude uncertainties can be 
taken to be approximately constant for all normal spiral 
galaxies. 

5.1. Comparisons with Alternative Analyses 

Our Archangel analysis procedures can be com- 
pared with alternative reductions of Spitzcr observa- 
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Fig. 7. — Variation of magnitude uncertainty as a function of the 
sky value in MJy sr — . Squares represent low surface brightness 
galaxies while triangles stand for irregular ones. Galaxies repre- 
sented by asterisks are galaxies with a very bright object close to 
them. Sky uncertainties can be 4 times higher than normal without 
resulting in abnormally high uncertainty in magnitudes. 

tions. Compariso ns with magnitudes found by the 
projects SINGS (iMu fioz-Mateos et~aTl [200l . S4G 
Sheth et al.l (12010J). pr ivate communication), and CHP 



Freedman et al.l |2011|) . private communication) show 



that our Cosmic Flows project is on the same magnitude 
scale as all these projects. In the case of CHP we 
give special attention to a comparison because our two 
programs, CHP and CFS, have the common ambition 
of measuring galaxy distances. As we go forward, we 
want to understand to what degree the alternative 
photometry analyses are interchangeable. A comparison 
is given between the two sources in Figure [9l There 
is a slight tendancy for CHP values to be brighter 
for the largest galaxies, with essentially no difference 
faintward of [3.6] = 12. The most likely explanation 
for a difference with the bright, large galaxies is small 
differences in the way sky values are set. The rms scatter 
in the differences (6 deviant points rejected) is ±0.052 
which, if attributed equally, implies an uncertainty in an 
individual measurement of ±0.037 mag for each source. 

Comparisons with other projects give comparable 
results. Typical zero point differences are ±0.01 and 
rms uncertainties arc ±0.04 — 0.05 mag. See a summary 
of comparisons with other major programs in Table [2j 
These results provide an estimate of the internal errors 
of alternate fitting procedures with the same data. We 
recall that our two measures of magnitude agree at the 
level of 0.01 with scatter ±0.02. 

In summary of errors, the dominant contributions are 
sky settings (0.04 mag), flux calibration (0.02 mag), and 
extinction (0.02 mag), leading to total uncertainties in 
magnitudes of ~ 0.05 mag. The great interest with the 
Cosmic Flows program is to use TFR to measure dis- 
tances to galaxies and the typical scatter in TFR is 0.4 
mag, 20% in distance. With errors on photometry af- 
ter corrections held to 0.05 mag the contribution to the 
distance error budget from photometry is minor. 

6. CONCLUSIONS 

Our photometric procedures for the semi-automated 
analysis of Spitzer IRAC channel 1 data at 3.6^m 
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Fig. 8. — Variation of magnitude uncertainty as a function of 
the galaxy: (top) semi-major axis, (middle) inclination, and (bot- 
tom) apparent area. Squares, triangles and asterisks represents 
low surface brightness, irregular, and possibly contaminated galax- 
ies respectively. There is no apparent correlation between magni- 
tude uncertainty and the radius, inclination, or apparent area of a 
galaxy. 

have been described. The galaxy surface photometry 
is ca rr ied out with the Ar c hange l software (|Schombertl 
I2007t iSchombert fc Smith! 12012(1 adapted for Spitzer 
data input. Already material is available for some 
3000 galaxies from the Spitzer Heritage Archive and 
our Cosmic Flows with Spitzer program will supply 
information for an additional 1274 galaxies. The 235 
galaxies analyzed in the course of this paper will be 
used in a subsequent paper for the calibration of the 
mid-infrared Tully-Fisher relation. 



[3.6] Photometry 



9 



1 


2 


3 


4 


5 


Program 


N 


Range 


lY^cfs Toother 


rms 


SINGS 


12 


8-10 


-0.02 


0.07 


Jarrctt 


5 


9-11 


-0.01 


0.03 


S4G 


29 


9-16 


-0.01 


0.04 


CHP 


171 


8-16 


0.01 


0.05 



TABLE 2 

Comparisons between CFS magnitudes and other 

MAGNITUDES: (1) PROGRAM NAME, (2) NUMBER OF GALAXIES 
COMPARED, (3) RANGE OF MAGNITUDES, (4) DIFFERENCE CFS 
MAGNITUDE — OTHER MAGNITUDE, AND (5) SCATTER. 
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Fig. 9. — Comparison of Archangel exponential disk extrapolated 
magnitudes [3.6] ex t with aysmptotic total magnitudes from the 
Carnegie Hubble Program [3.6]chp- The slight tilt of the thick 
line best fit, and the offset of < [3.6] el t — [3-6]chp >= 0-01 from 
the dotted line null value, have only 2 — 2.5cr significance. Both the 
tilt and the offset have been computed rejecting 6 deviant points. 



The final goal of our project is to measure distances, 
hence map peculiar velocities, across the local universe 
within 10, 000 km s _1 using the correlation between 
galaxy luminosities and their rotation rates. We have 
demonstrated the ability to use Spitzer Space Telescope 
mid-IR data to perform surface photometry with a 
relatively high accuracy. No correlation is found be- 
tween magnitude uncertainties and other important 
galaxy parameters such as inclination, apparent area, or 
semi-major axis. We conclude that, after all corrections, 
uncertainties on magnitudes arc of the order ±0.05 for 
the regular spiral galaxies at the heart of our project. 
These uncertainties are small compared with the overall 
scatter in the TFR. Low surface brightness galax- 
ies or very irregular ones require special attention but 
these classes of galaxies are not of principal interest to us. 

For our purposes, the advantages of mid-infrared pho- 
tometry from space include minimization of both galactic 
and internal obscuration issues, very low backgrounds, 
and source fluxes dominated by old stellar populations 
that are good representatives of the baryonic mass. The 



most outstanding advantage, though, is the integrity 
and consistency of the photometry in all quadrants of 
the sky. The Extragalactic Distance Database, EDD, 
contains HI profile information that provides useful 
line widths for over 11,000 galaxies. Ongoing Spitzer 
observations are providing the complementary photo- 
metric information required for a dense, detailed map 
of structure and motions in the near part of the Universe. 



Acknowledgements. Comparisons between alternative 
analysis procedures have been facilitated by data made 
available by Tom Jarrett, Mark Scibcrt representing 
CHP, the Carnegie Hubble Program, and Kartik Shcth, 
representing S4G, the Spitzer Survey of Stellar Structure 
in Galaxies. Thanks to James Shombert for the devel- 
opment and support of the Archangel photometry pack- 
age. The referee Michael Pohlen encouraged us to write 
the overview of Cosmic Flows which constitutes the Ap- 
pendix and was otherwise very helpful. NASA through 
the Spitzer Science Center provides support for CFS, 
Cosmic Flows with Spitzer, cycle 8 program 80072. In 
addition to the authors, CFS co-investigators are Wendy 
Freedman, Tom Jarrett, Barry Madore, Eric Persson, 
Mark Seibert, and Ed Shaya. RBT receives support for 
aspects of this program from the US National Science 
Foundation with award AST-0908846. 

Appendix: Cosmic Flows Program Overview 

Cosmic Flows may have as many arms as an octopus. 
At its core is a collaboration between Courtois and 
Tully to obtain accurate distances to galaxies. A major 
part of the program involves exploitation of the TFR. 
Activities in this regard began with the accumulation 
of HI profiles for the necessary kinematic information 
within the Cosmic Flows Large Program using the 
US National R a dio A s tronom y Green Bank Telescope 
(|Courtois et al.l 120091 l2011b[) and the accumulation 
of optical photometry for the necessary magnitude 
and inclination informat ion using the Unive rsity of 
Hawaii 2.24m telescope (jCourtois et al.l I2011al ). The 
present extension embarks on complementing the optical 
photometry with mid-infrared photometry The current 
paper describes analysis procedures developed within the 
core program but Cosmic Flows with Spitzer embraces 
the larger team identified in the acknowledgement. The 
next paper in this series will include most of the CFS 
team in a discussion of the TFR calibration with Spitzer 
photometry 

Near, intermediate, and far TFR s amples in the 
Cosmic Flows program were described bv lCourtois et all 
(|2011b[) . The 'near' sample is intended to achieve dense 
coverage of a volume extending to 3300 km/s with 
inclusion of all galaxies typed later than Sa that are 
brighter than Mk = — 21, inclined greater than 45°, and 
not obscured, disrupted, or confused. The 'intermediate' 
sample is drawn from flux and color limits applied 
to an Infrared Astro nomical Satellite redshift survey 
()Saunders et al.l 120001 ) . The flux limit at 60 (im is 0.6 
Jy, the color criterion to separate normal spirals from 
active nuclei is a ratio of 100 fim to 60 fim flux greater 
than one, there is a velocity cutoff at 6000 km/s, and 
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there is the same inclination restriction as with the 
near sample. By contrast, the 'far' sample is restricted 
to extreme edge-on systems drawn f r om Flat Galaxy 
catal ogues (jKarachentsev et al.l 119991 : iMitronova et alJ 
|2004| ). Candidates in the sample that lie at declinations 
accessible to Arecibo Telescope have velocities extending 
to 15,000 km/s. These are our well defined samples. 
In addition we derive distances to all other suitably 
observed galaxies. Generally the information for the 
additional systems comes from archives. In all, presently 
good data are available for about 7500 appropriate 
galaxies. 

A quite separate and active component of Cosmic 
Flows is a program with Hubble Space Telescope to 
obtain Tip of the Red Giant Branc h distances to 
nearby, spatially res o lved galaxies ([ Makaro v et alJ 
[200l 1miiUt~aIf [20071 IJacobs et all I2009D . Exquisite 
distances (5% accuracy) are available for approaching 
300 galaxies within ~ 10 Mpc. 

Distances for Cosmic Flows encompass measures 
by other methodologies discussed in the literature. 
Foremost among these are Cepheid Period-Luminosity 
Relation, Surface Brightness Fluctuation, Fundamental 
Plane, and Supernova la procedures. The diverse 
material is drawn tog ether in EDD, th e Extragalactic 
Distance Database^ (|Tullv et al.l [20091) . EDD goes 
beyond the compilation of catalogs relevant to extra- 
galactic distances to include redshift catalogs, that with 
various levels of completion describe the distribution of 
galaxies in the local universe, and group catalogs, that 
help identify entities where averaging over velocities or 
distances is reas onable. The first assembly of distances 
in this program (jTullv et al.ll200l has now been given 
the name Cosmicflows-1. A core team is now involved in 
the assembly of Cosmic flows-2 (jTullv &: Courtoid 120121 : 
ICourtois fc Tuhvll2012bfl . 

The holy grail of Cosmic Flows is the use of distances 
to determine peculiar velocities and, subsequently, mass 
fluctuations. Peculiar velocities are departures from the 
cosmic mean expansion and it is assumed that they 
arise due to density irregularities. Two regimes require 
separate attention. The high density environments in 
and around collapsed halos are at the extreme of non- 
linear dynamics. Within the collaboration we have de- 
veloped Numerical Action Methods that provide an op- 
timal description of the distribution of mass affecting 
gala xies on curved orb i ts on first approach to an attrac- 
tor (IShava et al.l [T995t IFeebles et al.l[200l [20I1 . The 
other extreme is the regime of linear dynamics. A pro- 
cedure we have used that is appropriate with redshift 
data sets of 10 5 or more obje cts is based on the action 
principle (jLavaux et al.ll20io() . However the methodol- 
ogy that most interests us starts with Wiener filtering 
of the peculiar velocity field resulting in descriptions of 
the density field in dependent of information provided 
by re dshift surveys (jZaroubi et al.l 119951 : ICourtois et alJ 
120121 1. The current density field can be mapped back 
to initial conditions that are then the starting point for 
constrained simulations that attempt to approximate the 



observed universe with a computer model dKlypin et all 
[2001 IGottloeber et alJl26Tol: ICourtois fc Tullvll2012aL 



5 http://edd.ifa.hawaii.edu 
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